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ABSTRACT: In a recent letter, Govaert et al. examined the relationship 
between strain hardening modulus and flow stress (Tfiow for five different 
glassy polymers. In each case, results for Gr at different strain rates or different 
temperatures were linearly related to the flow stress. They suggested that 
this linear relation was inconsistent with simulations. Data from previous 
publications and new results are presented to show that simulations also yield 
a linear relation between modulus and flow stress. Possible explanations for 
the change in the ratio of modulus to flow stress with temperature and strain 
rate are discussed. 
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INTRODUCTION 

The response of polymers to large strains plays a key role in determining their failure mode 
and mechanical performance. Strain hardening, an increase in stress with increasing strain, 
makes it more difficult to deform regions that have already yielded. This spreads deformation 
to new regions, preventing the localization of strain that would lead to brittle fracture. 

Many experimental measurements of strain hardening have been interpreted in terms of 
an entropic network model ^ based on rubber elasticity theories.- The idea is that the entan- 
glements between chains act like chemical crosslinks. The increase in stress is attributed to 
the decrease in entropy as the chain segments between these effective crosslinks are stretched. 
While this theory describes the functional form of stress-strain curves, the magnitude of the 
rate of stress increase, or hardening modulus Gr-, is hard to understand from this model.-i^ 
Assuming that the entanglement density is comparable to that in the melt, the hardening 
modulus should be comparable to the melt plateau modulus Gmeit near T^, and decrease 
with decreasing temperature. Instead Gr is typically two orders of magnitude larger than 
Gmeit and increases with decreasing T. 

In a series of recent papers,-i^i^ we have used molecular simulations to probe the origins 
of strain hardening. The microscopic behavior is not consistent with the entropic network 
model. Instead, strain hardening is directly related to the rate of plastic rearrangements 
needed to maintain chain connectivity. The scale of the hardening modulus Gr is thus set 
by the flow stress CTflow rather than entropic stresses. In both these simulations^>^!^ and 
experiments,-!^*^ the values of Gr and ciflow are of the same order of magnitude and both 
increase with decreasing temperature. Our simulations used a coarse-grained bead-spring 
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model,— but qualitatively similar behavior is observed with more realistic potentials.— A 
recent microscopic theory of glassy polymers also leads to similar scaling of plastic flow and 
hardening stresses.— 

Motivated by these studies, Govaert et al.— examined the hardening modulus and flow 
stress of five different glassy polymers over a range of rates and temperatures. They found 
a very interesting linear relationship between Gr and a flow 

Gr = Co + Cifffiow (1) 

for all the polymers studied. Unfortunately they took our statement that the yield stress 
sets the scale of the hardening modulus to mean that the ratio of hardening modulus to 
yield stress is strictly constant, and thus argued that the observation of a nonzero Gq seemed 
inconsistent with our results. This is not the case, but their work has inspired us to re- 
examine the relation between Gr and (Tfiow 

In this paper we compare previously published simulation results^"^ and new data to 
the experimental data shown by Govaert et al..— All simulations are consistent with Eq. [H 
and Co is not in general zero. The implications of Eq. [1] as the temperature T approaches 
the glass transition temperature Tg are discussed, as well as possible origins for Gq due to 
variations in thermal activation with strain. 

SIMULATION METHODS 

The simulations follow the methodology described in our previous papers.-i^>^ A generic bead- 
spring model ^'^ that describes the coarse-grained behavior of polymers is used. Each polymer 
contains = 350 spherical beads that interact with a truncated and shifted Lennard- Jones 
(LJ) potential Ulj = Auq [{a/ry^ — (a/r)^ — {a/vcY^ + (a/rc)^]. The binding energy uq and 
molecular diameter a are used to define our units. The unit of time r = ^yma'^/uo, where 
m is the bead mass. 

Adjacent beads along the chain are coupled with the finitely extensible nonlinear elastic 
(FENE) potential, UpENEir) = —kRlln{l — {r/Ro)"^), which prevents chain crossing and 
scission. The standard— values k = SOuq/o^ and Rq = 1.5a are employed, giving an equi- 
librium bond length Iq — 0.96a. The entanglement density is varied by adding a bending 
potential Ubend to change the chain stiffness. This potential, UbendiG) = fcfeend(l — cos 6*), 
where 9 is the angle between consecutive covalent bond vectors. Increasing /cfeend from to 
2.0mo changes the number of monomers per entanglement length N^, from about 70 to 20.— 

Cubic samples with Nch = 200 chains and periodic boundary conditions are equilibrated 
using the double-bridging-MD hybrid (DBH) algorithm.— A uniaxial compression is applied 
along the z axis while maintaining zero stress along the transverse (x, y) directions.— The 
stretch A = L^/L^, where Lz is the period along the z direction and L° is the initial value. A 
constant true strain rate e = A/A is applied. The stress a along the compressive axis is plotted 
against the Green-Lagrange strain g{\) = A^ — 1/A, since the slope of this curve corresponds 
to the hardening modulus in entropic models and experimental analysis: Gr = da/dg. 

The strain rates e used here range from |e| = 10~^/r to 10~^/r. Over this range we 
find no qualitative change in behavior, just a roughly logarithmic shift in stress with rate 
that is also seen in many experiments^ and is consistent with activated models^ for flow 
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stress. Since mappings of the bead-spring model to real systems give r in the picosecond to 
nanosecond range (e. g. 66ps for PE^), the slowest rates used here overlap with the highest 
rates accessible in experiments (10^ to 10^s~^).—^ 

In both simulations and experiments, the initial yield behavior is sensitive to the age 
and preparation of the glassy system.-i^^ Shear rejuvenates the system and the behavior for 
\g\ > 0.5 is fairly independent of past history. Experimental data typically show a large 
initial yield stress followed by strain softening. There is then a minimum stress before strain 
hardening sets in. Govaert et al. associated this minimum stress with the flow stress.— The 
value of \g\ at the minimum varies between about 0.5 and 1 with polymer and temperature, 
and the stress is fairly constant over this range of \g\. In our simulations, the glass has not 
generally been aged long enough to produce significant strain softening. There is an initial 
elastic increase in stress followed by a nearly constant plateau. We identified the stress at a 
point near the end of this plateau \g\ = 0.5 with the flow stress. This is close to the minimum 
in the stress for simulations that do show strain softening. 

Since the experiments and simulations deal with compression, the stress along the com- 
pressive axis is negative. Following Govaert et al. we will take Cfiow as a positive number 
equal to the magnitude of the stress. As in Ref.-, the value of Gr was obtained from a linear 
fit to 0" as a function of g from g = —0.5 to -3. 

RESULTS AND DISCUSSION 

Figure d] shows plots of a/ajiow for several temperatures and a fixed er = —3.16 ■ 10""^. Note 
that these normalized curves show an increase in slope with increasing T, implying that the 
ratio Gr/(Jfiow increases with increasing temperature. Using a value of |5'(A)| other than 0.5 
to define the flow stress would not change this trend. Govaert et al. found a similar increase 
for poly (ethylene terphthalate)-glycol (PETG) and a smaller increase for polystyrene (PS).— 
The trend with temperature had the opposite sign for polycarbonate (PC) and poly(methyl 
methacrylate) (PMMA). 

Given the observed change in Gr/c"flow, Govaert et al. examined the functional relation 
between the two quantities in more detail and discovered that results for all polymers could 
be fit to Eq. [H Figure [2] shows that our simulation results follow the same linear relation. 
Results for the temperature dependence at two values of chain stiffness are shown by solid 
symbols. The two have very different slopes, but similar positive offsets. We found that 
decreasing the strain rate by a factor of 30 only changed the offset by about 40%. As 
expected from entropic models, and our picture of hardening through plastic rearrangements, 
the slope is always steeper for more entangled polymers (bigger kbend)- Note that the unit 
of stress^ Uq/q^ is of order 50MPa so that the simulation results for Or and a flow are quite 
comparable to the experimental values (5 to 40MPa and 15 to 120MPa, respectively). 

The khend = 0.75mo data at fixed shear rate shown in Fig. [2] (circles) was included in 
Table 2 of Ref.^. This table also quoted values of Gr/cfiow, which systematically increased 
by 50% with increasing T. Given this, it is surprising that Govaert et. al. concluded that 
our simulations gave a constant Gr/cJfiow Note that the ratio between Gr and the entropic 
prediction changes by a factor of 50 for the same data. This is one reason for our conclusion 
that the scale of G^ is set by (Jfiow rather than the entropic stress. 
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Figure 1: (Color online) Compressive stress —a normalized by Ojiow as a function of —g{\) 
for systems with kbend = 0.75mo at ksTje = 0.01 (solid), 0.1 (dotted), 0.2 (dashed) and 0.3 
(dash-dotted). The strain rate e = —3.16 ■ 10~'^r~^. (Rescaled data from Fig. 3 of Ref. 5.) 
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Figure 2: (Color online) Hardening modulus as a function of flow stress for systems with 
kbend = 0.75mo (circles) and 1.5uo (squares). Filled symbols show results for fixed strain rate 
and four values of the temperature, ksT/uo = 0.01, 0.1, 0.2 and 0.3. Open symbols show 
results for kBT/uo = 0.2 and six values of strain rate |e|r = 10^^, 10^^, 3.16 • 10"^, 10"'', 
3.16- 10~^, and 10^'^. Both Gr and ctaow decrease with increasing temperature and decreasing 
strain rate. The straight lines are linear fits like those obtained in Ref.—. 
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The linear relation between Gr and (jfiow leads to interesting behavior as T increases to 
the glass transition temperature Tg where a flow —^0. In Fig. [21 the results for varying T 
have Co > 0, implying that hardening would still be observed in the absence of a flow stress. 
While a hardening modulus of order the melt plateau modulus would be expected in this 
limit, the values of Co in our simulations and for PETG and PS are substantially larger than 

One possible explanation of the large Cq was mentioned in the discussion of changes in 
Gr jo flow in Ref.-. Strain hardening is correlated with an increase in the size of regions 
that must be plastically deformed to shear the system while retaining chain connectivity. 
This is evidenced by an increase in the magnitude of non-affine deformations^ and in the 
number of Lennard Jones bonds broken-"^ as \g\ increases. The flow stress goes to zero 
when T is high enough to activate local segment-scale rearrangements at the given strain 
rate. This temperature will not in general be high enough to activate the larger, correlated 
rearrangements needed to maintain chain connectivity at higher strains, leading to a nonzero 
Gr as (Jficy^ goes to zero. 

Govaert et al. also found cases where fits to the temperature dependent hardening gave 
negative Cq. This implies that the hardening modulus vanishes and then becomes negative as 
T approaches Tg. This behavior seems counterintuitive, but the only case where the fit clearly 
indicates a negative Co is PC, and the data for it show a systematic curvature towards the 
origin. It would be interesting to extend the experimental data for all these systems toward 
Tg. One possibility is that the strain softening contribution to a is still important at the 
minimum in the stress that was used to define experimental values of Cfiow Indeed, the 
minimum occurs where the strain softening and straing hardening terms balance. The strain 
softening term in the experimental data varies in magnitude and extent with temperature, 
as does the value of \g\ at the stress minimum. 

The rate dependence of G/j is also of interest. In Ref.- we found that plots of a j a flow for 
different strain rates and fixed temperature collapsed fairly well onto a single curve, implying 
a small Cq. Govaert et al. found fairly good collapses in similar plots for PS and PMMA, 
but a clear change in normalized stress curves for PETG, PC and poly(phenylene ether) 
(PPE). All but PETG had a negative Cq. Inspired by the experiments, we have extended 
our results to strain rates that are thirty times lower than in Ref.-. As shown by the open 
symbols in Fig. [2l results for C^ vs. Ofiow fall on a straight line with a negative offset. Ref.- 
only included the four highest rates (and thus highest a flow)- While the ratio of Gr/<Jfiow 
changes relatively little for these four points, even these data are consistent with a negative 
Co. 

Rate dependence is often discussed in terms of Eyring-like models^i^S where the stress 
rises with rate as {ksT/V*) Ine. Here V* is called the activation volume and represents the 
derivative of the activation energy barrier with respect to stress. Both experiments-'- and 
simulations ^'^^ show a logarithmic dependence of stress on rate, with a strain dependent 
prefactor that increases (implying V* decreases) during strain hardening^^. This would be 
consistent with the complex rearrangements at large strains being harder to activate with a 
simple compressive stress. The more rapid rate dependence at high strains leads to a drop 
in strain hardening with decreasing rate, and thus is consistent with negative values of Cq. 
The only system considered by Govaert et al. that did not have a negative Co was PETG, 
which showed little change in either Gr or a flow- Note that a negative Co is not problematic 
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for the rate dependence since it is not generally possible to bring either Gr or a flow to zero 
because of the weak logarithmic dependence on rate. 

Govaert et al. state that the variation of Grjofiow in experiments represents evidence for 
a "so-called back-stress, an elastic contribution to the strain-hardening response." They did 
not explain this conclusion, but cite other observations that are frequently cited as evidence 
for a back-stress. For example, polymers exhibit a memory effect, returning to nearly their 
undeformed shape when heated above Tg. In Ref.- we showed that simulations reproduce 
this effect.'^ They also allow the magnitude of the stress driving this recovery to be evaluated 
directly. We found that the stress is entropic in origin and only of order the melt plateau 
modulus, rather than the much larger Gr- 

Govaert et al. also argue that the Bauschinger effect in oriented polymer glasses implies 
a back stress. There is an asymmetry in the stress needed to deform the system, with a 
larger stress required in the direction that would lead to an increase in orientation. This 
phenomenon is also consistent with our simulations and physical picture. When an oriented 
system is deformed in a way that reduces orientation, the constraints of chain connectivity 
can be satisfied in many ways. There is no need for the correlated rearrangements that 
occurred during orientation and the stress will be smaller. 

CONCLUSIONS 

Simulations were used to examine the relation between the flow stress and hardening modulus 
as temperature or strain rate was varied. Contrary to statements in Ref.—, the results are 
completely consistent with experiments.— Results for both simulation and experiments show 
a linear relation between G^ and a flow In both cases the slope G\ and intercept Cq depend 
on whether T or e is varied and Cq is more positive for varying temperature than varying 
rate. The observations also seem to be more easily understood in terms of a strain hardening 
model based on plastic deformation than an entropic model. 

A qualitative explanation of the trends with temperature and rate was given based on 
the previously observed correlation between increases in stress and increases in the rate of 
plastic rearrangement .-1^1^ Since the size and complexity of plastic rearrangements increase 
with strain, thermal activation may be more effective for the initial flow stress than G^. This 
would lead the flow stress to vanish at a lower temperature than G^, implying a positive 
Co for variations with T. In contrast, the value of Co would be negative for rate-dependent 
data if the simpler activations involved in a flow lead to a stronger rate dependence. The 
simulations are consistent with these predictions for Co and most polymers have Co < 
for rate-dependent data. Larger values of Co are found in experimental results for temper- 
ature dependence, although in some cases they appear to be negative. This may reflect the 
contribution of strain softening to the measured values of Ofiow- 

Many aspects of strain hardening result from the orientation of molecules that is produced 
by strain. Any mechanism that depends on this orientation will have memory effects and 
may exhibit a Bauschinger effect, but different models give very different predictions for the 
magnitude of these effects. In conventional entropic network models, orientation represents a 
loss of entropy that leads to a direct entropic stress. This entropic term is much smaller than 
the hardening modulus, but comparable to the stress driving the shape memory effect.- Our 
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simulations suggest that the influence of orientation on Gr is indirect. As the orientation 
increases, larger non-afline displacements and more breaking of van der Waals bonds are 
required to maintain chain connectivity. The growing constraints of chain connectivity on 
rearrangements are related to the decreasing entropy. However the magnitude of the stress 
is much higher than the plateau modulus because the energy dissipated by rearrangements 
is set by the flow stress rather than /c^T. 
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